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Abstract 

The momentum transfer dependence of the total cross section for 
elastic scattering of cold dark matter candidates, i.e. lightest su- 
permymmetric particle (LSP), with nuclei is examined. We find that 
even though the energy transfer is small (< WOKeV) the momentum 
transfer can be quite big for large mass of the LSP and heavy nuclei. 
The total cross section can in such instances be reduced by a factor 
of about five. 



There is ample evidence that about 90% of the matter in the universe 
is non-luminous and non-baryonic of unknown nature [1-3]. Furthermore, in 
order to accommodate large scale structure of the universe one is forced to 
assume the existence of two kinds of dark matter [3]. One kind is composed 
of particles which were relativistic at the time of the structure formation. 
This is called Hot Dark Matter (HDM). The other kind is composed of par- 
ticles which were non-relativistic at the time of stucture formation. These 
constitute the Cold Dark Matter (CDM) component of the universe. The 
COBE data [4] by examining the inisotropy on background radiation suggest 
that the ratio of CDM to HDM is 2:1. Since about 10% of the matter of the 
universe is known to be baryonic, we know that we have 60% CDM, 30% 
HDM and 10% baryonic matter. 

The most natural candidates for HDM are the neutrinos provided that 
they have a mass greater than leV/c 2 . The situation is less clear in the case 
of CDM. The most appealing possibility, linked closely with Supersymmetry 
(SUSY), is the LSP i.e. the Lightest Supersymmetric Particle. 

In recent years the phenomenological implications of Supersymmetry are 
being taken very seriously [5-7]. Pretty accurate predictions at low energies 
are now feasible in terms of few input parameters in the context of SUSY 
models without any commitment to specific gauge groups. More or less such 
predictions do not appear to depend on arbitrary choices of the relevant 
parameters or untested assumptions. 

In such theories derived from Supergravity the LSP is expected to be 
a neutral fermion with mass in the 10 — 100GeV/c 2 region travelling with 
non-relativistic velocities (f3 ~ 10~ 3 ) i.e. with energies in the KeV region. In 
the absence of R-parity violation this particle is absolutely stable. But, even 
in the presence of R-parity violation, it may live long enough to be a CDM 
candidate. 

The detection of the LSP, which is going to be denoted by xi, is extremely 
difficult, since this particle interacts with matter extremely weakly. One 
possibility is the detection of high energy neutrinos which are produced by 
pair annihilation in the sun where this particle is trapped i.e. via the reaction 

Xi + Xi^v + v (1) 
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The above reaction is possible since the LSP is a majorana particle, i.e. its 
own antiparticle (a la 7r°). Such high energy neutrinos can be detected via 
neutrino telescopes. 

The other possibility, to be examined in the present work, is the detection 
of the energy of the recoiling nucleus in the reaction 



This energy can be converted into phonon energy and detected by a tem- 
perature rise in cryostatic detector with sufficiently high Debye temperature 
[3,8,9]. The detector should be large enough to allow a sufficient number of 
counts but not too large to permit anticoincidence shielding to reduce back- 
ground. A compromise of about lKg is achieved. Another possibility is the 
use of superconducting granules suspended in a magnetic field. The heat 
produced will destroy the superconductor and one can detect the resulting 
magnetix flux. Again a target of about lKg is favored. 

There are many targets which can be employed. The most popular ones 



contain the nuclei \He, fF, ™Na, f Ca, 72 f 2 Ge, 7 £As, If I, IfXe, and 



It has recently been shown that process (2) can be described by a four 
fermion interaction [10-16] of the type [17] 



Xl + (A,Z)^ X i + (A Z) 



(2) 





G_f_ 

V2 



[</aXi7 A 7 5 Xi + JxiXi] 



(3) 



where 



Jx = Nlxify + f v n + {f A + /Ar 3 ) 7 5 }N 



(4) 



and 



J = N(fs + fs^)N 



(5) 
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where we have neglected the uninteresting pseudoscalar and tensor currents. 
Note that, due to the majorana nature of the LSP, Xil X Xi — (identically). 
The vector and axial vector form factors can arise out of Z-exchange and 
s-quark exchange [10-15] (s-quarks are the SUSY partners of quarks with 
spin zero). They have uncertainties in them (see ref. [15] for three choices 
in the allowed parameter space of ref. [5]). The transition from the quark 
to the nucleon level is pretty straightforward in this case. We will see later 
that, due to the majorana nature of the LSP, the contribution of the vector 
current, which can lead to a coherent effect of all nucleons, is suppressed [10- 
15]. Thus, the axial current, especially in the case of light and intermediate 
mass nuclei, cannot be ignored. The scalar form factors arise out of the 
Higgs exchange or via S-quark exchange when there is mixing between s- 
quarks q L and q R [10-12] (the partners of the left-handed and right-handed 
quarks). They have two types of uncertainties in them [18]. One, which is 
the most important, at the quark level due to the uncertainties in the Higgs 
sector. The other in going from the quark to the nucleon level [16-17]. Such 
couplings are proportional to the quark masses and hence sensitive to the 
small admixtures of qq (q other than u and d) present in the nucleon. Again 
values of fg and f s in the allowed SUSY parameter space can be found in 
ref. [15]. 

The invariant amplitude in the case of non relativistic LSP takes the form 

[15] 



|m| 2 = ^- mi 2 +P *" P V o| 2 +|J| 2 + |J| 2 

~ /? 2 |Jo| 2 + |J| 2 + l^| 2 (6) 

where | Jq\ and |J| indicate the matrix elements of the time component and 
space component of the current | J\\ of eq. (4) and J the matrix element of 
the scalar current J of eq. (5). Notice that |J | 2 is multiplied by f3 2 (the 
suppression due to the majorana nature of LSP mentioned above). It is 
straightforward to show that 



\M 2 = A*\F(e)\ 2 (ft - ti^-J (7) 
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J 2 = A 2 |F(q 2 )| 2 (f s -fl^) 2 (8) 
|J| 2 = 1 < J iW [/^o(q) + /A^i(q)] ||4 > I 2 (9) 

with 



A A 

n (q) = £ "(/)' ; " x '- ni(q) = E ^0>sO>" iq ^ (io) 

where cr(J), r 3 (j), Xj are the spin, third component of isospin (r 3 \p >= \p >) 
and cordinate of the j-th nucleon and q is the momentum transferred to the 
nucleus. 

The differential cross section in the laboratory frame takes the form [15] 



da a 0/ mi 2 1 fffl2|f|2[1 2?? + 1 2 2 2 

^ = (TT^ {/3 1 Jo1 [1 " JT+W ] + 1 1 + 11 } ( } 

where r] = m\jm v A ( m p = proton mass), (3 = v/c (v is the velocity of LSP), 
mi is the mass of LSP, £ = • q > (forward scattering) and 

a = -L(G F m p f ~ 0.77 x l(T 38 cm 2 (12) 

| Jo| 2 , |J| 2 and | J| 2 are given by eqs. (7)-(9). The momentum transfer q is 
given by 

|q| = QoC, Qo = fi 2 ^- (13) 

Some values of go (forward momentum transfer) for some characteristic 
values of m\ and representative nuclear systems (light, intermediate and 
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heavy) are given in table 1. It is clear that the momentum transfer can be 
stable for large mi and heavy nuclei. 



The total cross section can be cast in the form 



+ {ft - fl^f Wo) - Y^TW^ f y - fv^^o)] 

+ (/^ (0)) 2 /oo(g 2 ) " 2/^ (0)fi 1 (0)/ 01 (g 2 ) 

+ (/A^i(0)) 2 /ii(g 2 )} (14) 

where 

/p(% 2 ) = 2(p + l) fe +2p \F{qie)\ 2 d^ p = 0,l (15) 
j o 

n p (q) = (2J, + < Jj| |O p (q) 1 1 Jj >, p = 0, 1 (16) 
(see eq. (10) for the definition of Vt p ) and 

V(9o)-2y o ^ -yyoT "n^oT P,P -0,1 (17) 

In a previous paper [16] we have shown that the nuclear form factor can 
be adequately described within the harmonic oscillator model as follows 

F(q 2 ) = [ jHqb, Z) + -<S>(qb, N) ] e"^/ 4 (18) 
where $ is a polynomial of the form [18] 

N ma x(a) 

<&(?&,«) = £ 6<?\qb) 2 \ a = Z,N (19) 

A=0 
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N ma x{Z) and N max (N) depend on the major harmonic oscillator shell occu- 
pied by protons and neutrons [16], respectively. The integral I p {q^) can be 
written as 

I p {ql) -> I p {u) = \F(2x/b 2 )\ 2 dx, (20) 

J 

where 



u = q 2 b 2 /2, b = 1.0A 1/3 fm (21) 
With the use of eqs. (18), (19) we obtain 

I p {u) = ±{ Z 2 lf z [u) + («) + m^ N [u)} (22) 

where 



(p) _ W ^ flf) 00) 2 A +^(A + z/ + p)! ri _ u A t^ +p m k 

A=o v=0 up a K=0 is,. 

(23) 

with a,p = N,Z 

The coefficients 8^ for light and medium nuclei have been computed in 
ref. [16]. In table 2 we present them by including in addition those for heavy 
nuclei. The integrals I p {u) for three typical nuclei ( io^a, HGe and f^Pb ) 
are presented in fig. 1 as a function of mi. We see that for light nuclei the 
modification of the cross section by the inclusion of the form factor is small. 
For heavy nuclei and massive mi the form factor has a dramatic effect on the 
cross section and may decrease it by a factor of about 5. The integral h(u) 
is even more suppressed but it is not very important. 

The spin matrix elements depend on the details of the structure of the 
nucleus considered. So is the spin form factor. The spin matrix element, 
since it does not show coherence, is expected to be more important in the 
case of odd light and intermediate nuclei. In the present work we will examine 
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the q 2 dependence of the spin matrix element in the cases of l 2 7 Pb and l 9 F 
whose structure is believed to be simple. 



To a good approximation [15,17] the ground state of the g 2 Pb nucleus 
1/2 neutron hole in the g 2 s 



can be described as a 2si/ 2 neutron hole in the g 2 8 -P& closed shell. One then 



finds 

flo(q) = (l/>/3)F 28 (q 2 ), ^(q) = -(l/VS)F 2s (q 2 ) (24) 



and 



'oo — J oi 



hi =2 f ms{q 2 )? (25) 
Jo 



Even though the probability of finding a pure 2si/ 2 neutron hole in the \ 
ground state of 2 % 2 Pb is greater than 95%, the ground state magnetic moment 
is quenced due to the 1 + p-h excitation involving the spin orbit partners. 
Hence we expect a similar suppression of the isovector spin matrix elements. 
Thus we write 



|(l/2)' > gs = C \(2s 1 / 2 )~ 1 > [l + C 1 |[0i 11/2 (n)(0i 13 /2)~ 1 (n)]l + > 

+ C 2 \[^ /2 { P ){Qh 11/2 )-\ P )\l + >+...] (26) 
Retaining terms which are most linear in the coefficients Ci, C 2 we obtain 



fio(q) = C 2 {F 2s (g 2 )/v / 3-8[(7/13) 1 / 2 C 1 F 0j (g 2 ) + (5/ ll) 1 / 2 C 2 F 0h (q 2 )] } 

(27) 



fi!(q) = Cl{F 2s {q 2 )/^-W/^f' 2 CA{q 2 ) - {b / ll) 1 ' 2 C 2 F Qh {q 2 )]} 

(28) 

where 



Nmax 

F n M 2 )=e^ b2/i £ 7f } (^) 2A (29) 

A=0 
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The coefficients 7^ are given in table 3. 

The coefficients Co, C\ and C2 were obtained by diagonalizing the Kuo- 
Brown G-matrix [18,19] in a model space of 2h-lp configurations. Thus we 
find 



C = 0.973350, d = 0.005295, C 2 = -0.006984 
We also find 

n (0) = -(1/^(0.95659), (small retardation) (30) 

fti(0) = -(1/^(0.83296), [sizable retardation) (31) 

The amount of retardation of the total matrix element depends on the 
values of f\ and f\. Using eqs. (25) and (26) we can evaluate the integrals 
I 00 , ^01 an d ho- The results are presented in fig. 2. We see that for a heavy 
nucleus and high LSP mass the momentum transfer dependence of the spin 
matrix elements cannot be ignored. 

In the second example we examine the spin matrix elements of the light 
nucleus l 9 F. Assuming that the ground state wave function is a pure SU (3) 
state with the largest symmetry i.e / = [3], (A/i) = (60), we obtain [20,21] 
the expression 

§g = + §W> (32) 

There is only isovector component contribution (the isoscalar matrix el- 
ement vanishes). The results for the I\i{u) integral are shown in fig. 3. We 
see that the effect of the nuclear form factor on the cross section for light 
nuclei is insignificant. 

In the present paper we have examined the momentum transfer depen- 
dence of the nuclear matrix elements entering the elastic scattering of cold 
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dark matter candidates (LSP) with nuclei. We have found that such a mo- 
mentum transfer dependence is very pronounced for heavy nuclear targets 
and mass of the LSP in the 100 GeV region. 
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Figure Captions 



Fig. 1. The integral Io(u), which describes the main coherent contribution 
to the total cross section as a function of the LSP mass (mi), for three typical 
nuclei: f Ca, f 2 Ge and jfPb. 

Fig. 2. The integral I\(u), entering the total coherent cross section as 
a function of the LSP mass (mi), for three typical nuclei: 2oCa, l%Ge and 
IfPb. For its definition see eqs (11) and (15) of the text. 

Fig. 3. The integral In, associated with the spin isovector - isovector 
matrix elements for ™ Pb and l 9 F as a function of the LSP mass (mi). The 
other two intergals Too and 7 i are almost identical and are not shown. 
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Table 1 : The quantity go (forward momentum transfer) in units of fm 
for three values of mi and three typical nuclei. 



Nucleus 


mi = 30GeV 


mi = lOOGeV 


mi = 150GeV 


20 Ca 


.174 


.290 


.321 


llGe 


.215 


.425 


.494 


™Pb 


.267 


.685 


.885 



Table 2. The coefficients 9\ determining the proton and neutron form factors 
for all closed (sub)shell nuclei. In a harmonic oscillator basis they are rational 
numbers. The coeficients for A = are equal to Z (or N). 



nlj-level 


A = 


A = 1 


A = 2 


A = 3 


A = 4 


A = 5 A = 6 


OSi/2 
0^3/2 


2 
6 


-2/3 










Opi/2 


8 


-1 










04/2 


14 


-3 


1/10 










16 


-11/3 


11/60 








04/2 


20 


-5 


1/4 








0/7/2 


28 


-9 


13/20 


-1/105 






1P3/2 


32 


-11 


61/60 


-11/420 






0/5/2 


38 


-14 


79/60 


-1/30 






lPl/2 


40 


-15 


3/2 


-1/24 






0^9/2 


50 


-65/3 


5/2 


-5/56 


1/1512 




0^7/2 


58 


-27 


33/10 


-107/840 


1/840 




14/2 


64 


-31 


17/4 


-173/840 


1/336 




14/2 


68 


-101/3 


293/60 


-31/120 


1/240 




2si/ 2 


70 


-35 


21/4 


-7/24 


1/192 




0h 11/2 


82 


-45 


29/4 


-73/168 


37/4032 


-1/27720 


0h 9 /2 


92 


-160/3 


107/12 


-31/56 


151/12096 


-1/15120 


1/7/2 


100 


-60 


217/20 


-653/840 


449/20160 


-1/5040 


1/5/2 


106 


-65 


123/10 


-397/420 


199/6720 


-1/3360 


2^3/2 


110 


-205/3 


403/30 


-153/140 


253/6720 


-1/2240 


2pi/ 2 


112 


-70 


14 


-7/6 


1/24 


-1/1920 


0^13/2 


126 


-84 


35/2 


-3/2 


1/18 


-49/63360 1/617760 
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Table 3. The coefficients 7^ , entering the polynomial describing the form 
factor (see eq. (29)) of a single particle harmonic oscillator wave function up 
to 6hu, i.e. throughout the periodic table. 



n I 


A = 


A = 1 


A = 2 


A = 3 


A = 4 


A = 5 


A = 6 





1 














1 


1 


-1/6 












1 


1 


-1/3 


1/24 










2 


1 


-1/3 


1/60 










1 1 


1 


-1/2 


11/120 


-1/240 








3 


1 


-1/2 


1/20 


-1/840 








2 


1 


-2/3 


11/60 


-1/60 


1/1920 






1 2 


1 


-2/3 


19/120 


-11/840 


1/3360 






4 


1 


-2/3 


1/10 


-1/210 


1/15120 






2 1 


1 


-5/6 


17/60 


-31/840 


9/4480 


-1/26880 




1 3 


1 


-5/6 


29/120 


-47/1680 


37/30240 


-1/60480 




5 


1 


-5/6 


1/6 


-1/84 


1/3024 


-1/332640 




3 


1 


-1 


17/40 


-31/420 


27/4480 


-1/4480 


1/322560 


2 2 


1 


-1 


2/5 


-1/15 


41/8064 


-1/5760 


1/483840 


1 4 


1 


-1 


41/120 


-1/20 


1/315 


-1/11880 


1/1330560 


6 


1 


-1 


1/4 


-1/42 


1/1008 


-1/55440 


1/8648640 
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